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AMtlfl£l 

A  series  of  numerical  experiments  to  predict 
the  evolution  of  the  Gulf  Stream  have  been 
performed.  The  forecast  system  uses  a  primi¬ 
tive  equation  model  of  the  north  western  Atlan¬ 
tic  Ocean  and  assimilation  schemes  which  em¬ 
ploy  both  a  feature  model  and  statirica! 
correlations  derived  from  the  regional  climutol- 
ogy  of  in  situ  data  and  long  time-base  numeri¬ 
cal  simulations.  The  evaluation  criterion  is  the 
mean  absolute  distance  between  forecast  loca¬ 
tions  of  the  Gulf  Stream  front  and  actual  loca¬ 
tions  as  verified  from  extensive  satellite  and  in 
situ  data.  Eight  one-week  and  five  two-week 
evaluation  intervals  during  1986-1988  were  se¬ 
lected  to  represent  a  variety  of  both  active  and 
inactive  Gulf  Stream  regimes.  To  insure  ob¬ 
jectivity.  hindcasting  was  disallowed.  The  re¬ 
sulting  forecasts  were  significantly  better  than 
persistence  at  both  one-  and  two-week  inter¬ 
vals.  This  study  indicates  the  feasibility  of 
Gulf  Stream  forecasting  using  assimilation 
schemes  which  provide  adequate  deep  infor¬ 
mation  and  numerical  models  which  are  de¬ 
signed  to  be  consistent  with  available  data. 


1  Introduction 

The  western  boundary  current  regions  of  the 
oceans  represent  domains  of  high  variability 
and  considerable  eddy  activity.  They  are  not 
adequately  sampled  by  in  situ  measurements  or 
by  satellite  altimetry  ( 10,13]  to  pemiit  accurate, 
instantaneous  estimates  of  sea  surface  topogra¬ 
phy.  They  are  also  regions  of  significant  stra¬ 
tegic  importance.  The  development  of  numer¬ 
ical  models  of  the  North  Atlantic  1 1 5]  capable 
of  reproducing  the  measured  variabilities, 
large  scale  circulations,  and  eddy  activity  has 
given  us  the  confidence  to  proceed  with  the 
construction  of  a  system  which  will  use  such 
models  to  dynamically  interpolate  the  asynop- 
tic  measurements  available  and  thereby  provide 
both  improved  nowcast  and  forecast  capabili¬ 
ties. 

A  system  for  constructing  initial  state  estimates 
and  performing  forecasts  was  constructed  for 
the  region  of  the  North  Atlantic  by  the  Data  As¬ 
similation  Research  and  Transition  (DART) 


team  at  the  U.  S.  Navy's  Naval  Oceanographic 
and  Atmospheric  Research  Laboratory. 

In  order  to  provide  a  'level  playing  field'  in  the 
arena  of  ocean  forecasting,  a  set  of  reference 
data  (i.e..  Gulf  Stream  frontal  locations)  was 
created  for  a  subset  of  the  Gulf  Stream  region. 
Before  any  forecasts  had  been  done,  data  over 
the  preceding  three  years  was  examined  to  at¬ 
tempt  to  define  several  periods  for  w  hich  accu¬ 
rate  positions  of  the  Gulf  Stream  axis  and  ed¬ 
dies  could  be  constructed  for  at  least  two 
weeks.  Clear  satellite  IR  imagery  was  re¬ 
quired,  and  in  most  cases  XBTs  and  GEOSAT 
altimetry  [3]  was  used  to  refine  the  locations  of 
the  features.  Four  time  periods  were  chosen  in 
1986  through  1988  which  permitted  eight  1- 
week  foiucast  and  verification  experiments  and 
five  2-week  experiments. 

The  primary  evolution  criterion  was  the  mean 
absolute  distance  between  forecast  locations  of 
the  Gulf  Stream  front  and  the  actual  (or  verifi¬ 
cation)  locations.  Forecast  error  was  comput¬ 
ed  as  the  average  absolute  offset  between  the 
forecasted  position  of  the  axis  with  the  position 
given  in  the  verification  state.  Persistence,  the 
assumption  of  no  motion  over  the  forecast  in¬ 
terval,  was  used  as  a  comparison  in  judging 
forecast  skill.  For  a  model  to  have  any  signifi¬ 
cant  skill  in  forecasting,  the  error  in  its  forecast 
(the  forecast  error)  must  be  less  than  the  error 
computed  by  simply  using  the  initial  state  as 
the  forecast  (the  persistence  error). 

To  examine  the  degree  to  which  this  small 
number  of  states  was  representative  of  the  Gulf 
Stream  as  a  whole,  persistence  errors  were 
computed  both  from  these  states  and  from  a 
year  of  the  operationally  produced  front  loca¬ 
tions.  The  distributions  were  compared  to  ver¬ 
ify  that  the  cases  chosen  were  representative. 

The  domain  was  divided  into  three  subregions: 
a  Western  Region,  extending  from  73°W  to 
66°W  longitude,  a  Central  Region,  extending 
from  66°W  to  59°W  longitude,  and  an  Eastern 
Region,  extending  from  59°W  to  53°W  longi¬ 
tude.  The  average  absolute  offset  error  of  each 
of  these  subregions  was  computed,  as  well  as 
the  overall  region,  extending  from  73°W  to 
53°W  longitude. 
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Figure  2  presents  the  results  front  one  of  the  2- 
wcek  forecast  experiments.  Measuring  the  er¬ 
rors  in  the  overall  region  from  73°W  to  53°W 
longitude,  persistence  represents  an  error  of  50 
km  at  1  -week  and  65  km  at  2-wceks.  By  com¬ 
parison,  the  forecast  provided  errors  of  only 
26  km  and  39  km  at  1-  and  2-weeks  respec¬ 
tively.  When  all  eight  1-week  and  five  2- week 
cases  are  examined,  it  was  found  that  the  sys¬ 
tem  provided  forecasts  which  were  better  than 
persistence  by  approximately  6  km  at  1-week 
and  10  km  at  2-weeks.  Standard  statistical 
techniques  were  used  to  verify  the  significance 
of  these  results. 


3  DART  iaJad  Nowcasl/Forecast 
System 

In  this  section,  we  provide  a  detailed  descrip¬ 
tion  of  the  DART  Gulf  Stream  run  stream.  The 
procedure  used  to  initialize  and  run  the  DART 
forecasts  begins  with  the  subjective  preparation 
of  an  initial  Gulf  Stream  frontal  location  map. 
This  manually  prepared  map  blends  frontal  lo¬ 
cation  information  contained  in  satellite  IR  im¬ 
agery,  satellite  altimetry,  and  any  available  BTs 
into  a  continuous  depiction  of  the  surface  fron¬ 
tal  location.  The  resulting  continuous  depic¬ 
tion  of  North  Wall  location  is  then  run  through 
a  feature  model  algorithm  in  the  OTIS  comput¬ 
er  program  which  provides  an  initial  state  esti¬ 
mate  of  the  dynamic  height,  which  is  then  con¬ 
verted  to  an  upper  layer  pressure  anomaly  (pi) 
in  the  NOARL  primitive  equation  circulation 
model.  Scaling  between  OTIS  dynamic  height 
and  model  pi  is  necessary  to  correctly  repre¬ 
sent  the  transport  in  the  model's  thick  upper 
layer.  This  initial  pi  field  is  instantaneously 
converted  in  an  accompanying  lower  layer 
pressure  field  (p2)  using  a  statistical  inference 
technique  based  upon  the  circulation  model's 
climatology.  Together,  pi  and  p2  are  then 
used  to  provide  for  geostrophic  initialization  of 
the  circulation  model.  For  a  brief  interval  im¬ 
mediately  following  these  "cold  start"  initializa¬ 
tions,  a  filter  is  applied  in  the  circulation  model 
run  to  nearly  eliminate  gravity  waves.  Finally, 
the  pl=0  contour  in  the  model's  forecast  state 
is  used  to  define  the  forecast  frontal  location 
for  direct  comparison  with  independent  verifi¬ 
cation  frontal  location  maps. 

Much  of  the  success  in  these  "cold  start" 
DART  forecasts  is  a  result  of  the  several  tech¬ 
niques  used  to  directly  (and  instantaneously) 
transfer  upper  layer  information  into  the  lower 
layer  of  the  circulation  model.  These  tech¬ 
niques  are  described  chronologically  in  the  fol¬ 
lowing  sections. 


3.1  OTIS  Feature  Model 

Maps  of  surface  topography  used  to  initialized 
the  circulation  model  and  to  verify  the  forecasts 
were  prepared  using  the  regional  OTIS,  devel¬ 


oped  primarily  at  the  Fleet  Numerical  Oceanog¬ 
raphy  Center  (FNOC)  [5]  with  contributions 
from  the  Naval  Oceanographic  Office  (NOO) 
and  the  Naval  Oceanographic  and  Atmospheric 
Research  Laboratory  [  1 ,2)  OTIS  is  a  data  qual¬ 
ity  control  and  interpolation  system  which 
combines  climatology,  maps  of  front  and  eddy 
boundaries,  MCSSTs  (satellite  IR  multi¬ 
channel  sea  surface  temperatures),  and  meas¬ 
ured  temperature  profiles  to  form  gridded 
three-dimensional  synoptic  thermal  analyses 
for  selected  ocean  regions.  A  reduced  set  of 
the  OTIS  system  capabilities  was  used  in  this 
study:  the  surface  topography  maps  were  pro¬ 
duced  using  maps  of  front  and  eddy  positions 
as  the  only  data  source.  The  OTIS  software 
interprets  these  maps  and  applies  module  f~r 
the  Gulf  Stream  front  and  eddies  to  form  a 
gridded  three  dimensional  field  of  temperature. 
Relative  dynamic  height  at  the  surface  is  then 
computed  directly  from  the  grid  of  temperature 
profiles  using  relationships  derived  from  anal¬ 
ysis  of  regional  historical  temperature  and  sa¬ 
linity  data  sets. 

All  data  for  each  analysis  date  were  combined 
onto  a  single  map  consisting  of  frontal  path 
segments  from  IR,  locations  of  front  and  ring 
crossings  from  altimetry,  and  AXBT  locations 
coded  according  to  water  type.  An  unbroken 
frontal  path,  from  the  western  to  the  eastern 
boundaries  of  the  model  donviin,  was  drawn 
by  hand  through  the  composite  data  set  and 
then  digitized,  and  ring  radii  and  center  loca¬ 
tions  were  extracted. 

OTIS  uses  parametric  models  of  the  Gulf 
Stream  front,  eddies,  and  the  ambient  back¬ 
ground  which  have  been  developed  from  a 
combination  of  historical  observations,  simple 
dynamical  models,  and  information  obtained 
from  published  studies. 

The  OTIS  system  assimilates  synthetic  obser¬ 
vations,  true  observations  taken  at  irregular 
times  and  positions,  and  climatology  to  form 
synoptic  maps  using  optimum  interpolation 
|7,4)  The  synthetic  observations  provide  a 
high-spatial  resolution  data  set  within  and  near 
fronts  and  eddies  where  observational  data  is 
often  too  sparse  to  resolve  these  features. 
Also,  the  synthetic  profiles  provide  subsurface 
information,  whereas  most  measurements  are 
made  only  at  the  surface  from  satellites.  This 
study  uses  only  synthetic  profiles  to  construct 
the  final  field  of  temperature. 

The  optimum  interpolation  results  in  a  three- 
dimensional  grid  of  temperature  covering  the 
domain  of  the  grid  of  the  circulation  model. 
Since  salinity  is  not  available  from  this  analy¬ 
sis,  dynamic  heights  at  the  surface  are  comput¬ 
ed  from  the  same  relationships  between  dy¬ 
namic  height  and  temperature  used  in  modcV.ng 
the  structure  of  rings.  The  root-mcan-square 
error  in  dynamic  height  computed  by  this 
method  is  about  0.065  dynamic  meters. 


—  4fi  — 


3-2  Circulation  Model 

The  circulation  model  used  in  this  study  is  do¬ 
cumented  in  Thompson  and  Hurlburt  [14], 
Murlburt  and  Thompson  (9|  and  Thompson 
and  Schmitz  115).  It  is  an  n-layer,  primitive 
equation  model  covering  the  region  from  78°W 
to  45°W  longitude  and  3()°N  to  45°N  latitude.  It 
includes  large-amplitude  bottom  topography. 
The  model  domain  was  chosen  so  that  the  vari¬ 
ability  in  the  location  of  the  Gulf  Stream  en¬ 
trance  into  the  domain  would  be  small.  Figure 
1  shows  a  plot  of  weekly  axis  locations  taken 
from  a  year  of  NHOC  boguses.  The  outlined 
box  in  the  plot  represents  the  circulate  model 
•  domain.  The  relatively  small  variability  in  the 
position  of  the  Gulf  Stream  axis  where  it  enters 
the  model  domain  permits  the  inflow  to  be 
specified  at  a  fixed  location  south  of  Cape  llat- 
teras.  The  version  of  the  model  used  in  this 
evaluation  included  two  layers,  with  a  deep 
western  boundary  current  (15)  supplied  by  an 
inflow  port  in  the  north  eastern  part  of  the  low¬ 
er  layer.  The  model  used  in  these  experiments 
is  on  a  spherical  grid  with  a  resolution  of  1/6 
degree  in  longitude  and  1/8  degree  in  latitude, 
which  represents  a  spatial  sampling  of  approxi¬ 
mately  14  km  in  each  direction  at  the  center  of 
the  grid. 

Since  layer  thickness  is  included  among  the 
model  variables,  fluctuations  of  the  pycnocline 
can  be  modelled  by  changes  in  the  depth  of  the 
interface  between  the  upper  and  lower  layers. 
This  permits  a  more  efficient  representation  of 
the  dominant  dynamical  modes  in  the  domain 
than  is  possible  with  a  model  which  uses  fixed 
thickness  levels.  This  was  deemed  of  particu¬ 
lar  importance  in  these  experiments,  due  both 
to  the  manner  in  which  we  initialize  the  lower 
layer  (described  in  detail  in  the  next  section) 
and  the  number  of  experiments  contemplated. 


3.3.  Statistical  Jn.f?reDC£.-af_Sab,lh.canq- 
cline  .Information 

Information  on  the  subthcrmocline  has  been 
shown  to  be  extremely  valuable  in  forecasts 
based  on  numerical  simulations  of  the  Gulf  of 
Mexico  (8,1 1 1  and  the  Gulf  Stream  (6,12). 

Ixing  model  simulations  are  used  to  derive  sta¬ 
tistical  relationships  between  the  subthermo- 
cline  pressure  (p2)  at  any  given  grid  point  in 
the  model  and  the  surface  pressure  (pi)  at  an 
array  of  grid  points.  Parameters  which  control 
this  derivation  are  chosen  to  maximize  the  skill 
in  estimating  the  lower  layer  pressure  in  an  in¬ 
dependent  dataset.  That  is,  coefficients  are  de¬ 
rived  from  one  run  of  the  model  and  are  used 
to  estimate  the  lower  layer  in  an  independent 
run. 

Figure  3  shows  an  example  of  using  these  co¬ 
efficients  to  estimate  the  lower  layer  pressure 
for  the  Gulf  Stream.  Note  that  while  the  pat- 
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tern  correlation  between  pi  and  p2  on  this 
model  day  is  only  .26,  the  correlation  exceeds 
.9  between  p2  and  the  estimate  of  p2  computed 
from  pi  by  the  statistical  inference  technique. 
The  coefficients  derived  from  these  lengthy 
model  simulations  are  applied  t<>  the  surface 
height  fields  produced  by  the  thermal  analysis 
to  provide  an  estimate  of  the  lowei  layer  pres¬ 
sure  field,  and  thus  the  pycnocline  depth 
anomaly  for  each  of  the  forecast  dates.  The 
absolute  accuracy  of  this  estimate  of  the  depth 
of  the  pycnocline  has  not  been  quantified,  but 
for  the  purposes  of  initializing  the  circulation 
model,  it  represents  lower  layer  information 
which  is  dynamically  consistent  with  the  upper 
layer  information  provided  by  OTIS. 

To  calculate  the  contribution  to  the  forecast 
skill  of  including  this  deep  pressure  informa¬ 
tion,  forecasts  were  made  using  the  present 
reference  datasets  and  using  three  alternatives 
for  defining  the  lower  layer  pressure  field  in 
the  initial  states.  Forecasts  made  using  either  a 
reduced  gravity  initialization  or  using  the  mod¬ 
el's  climatological  lower  layer  pressure  field 
resulted  in  a  larger  error  than  those  made  with 
the  statistical  inference  initialization,  generally 
by  about  5  km. 


3.4  Geostrophic  Velocity  Initialization 

The  remaining  model  variables,  the  u  and  v 
components  of  velocity  for  each  of  the  two 
layers  are  computed  geostrophically.  For  ex¬ 
ample,  in  the  upper  layer: 

k  x  fvt  =  -gVq 

where /  is  the  Coriolis  parameter  {f  =  2co  sin  6, 
where  (o  is  the  angular  velocity  of  the  earth's 
rotation  and  9  is  the  latitude),  1c  is  a  unit  verti¬ 
cal  vector,  is  the  geostrophic  component  of 

the  current  and  T|  is  the  free  surface  anomaly. 

3.5  .gravity  .IVaveJEiltec 

Despite  the  statistical  inference  of  the  lower 
layer  and  the  geostrophic  velocity  initialization, 
there  will  inevitably  still  be  some  dynamic  im¬ 
balances  in  the  initial  state.  One  advantage  of 
the  primitive  equation  model  approach  is  that 
such  imbalances  will  be  converted  to  short- 
period  gravity  waves,  which  can  easily  be  re¬ 
moved  by  selective  filtering.  For  the  particular 
domain  of  this  model,  the  dominant  gravity 
wave  period  is  approximately  8  to  9  hours. 
These  waves  are  attacked  in  the  NOARL  circu¬ 
lation  model  by  a  time-domain  running  average 
with  a  length  of  8.5  hours  which  is  applied 
once  at  +12  hours  into  the  forecast  and  again  at 
+24  hours,  after  which  no  further  filter  is  done 
during  the  remaining  two  weeks. 


4  ..Sian  mar  £ 

A  series  of  numerical  experiment  to  predici 
the  evolution  of  the  Gulf  Stream  have  been 
performed.  The  nowcast/forecast  system  de¬ 
veloped  by  the  DART  team  at  NOARL  is  the 
first  such  system  to  show  significant  skill  in 
forecasting  the  evolution  of  the  Gulf  Stream 
frontal  axis,  providing  estimates  of  the  front 
that  are  (on  average)  6  km  better  than  persis¬ 
tence  at  1-week  and  10  km  better  at  2- weeks. 
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Figure  1  —  One  year  of  weekly  Gulf 
Stream  axis  positions.  Inner  box  rep¬ 
resents  the  domain  of  the  forecast  sys¬ 
tem. 
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